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Abstract 

We discuss status of the singularity problem in General Relativity and argue that the requirement 
that a physical solution must be completely free of singularities may be too strong. As an example, 
we consider properties of the integrable singularities and show that they represent light horizons 
separating T-regions of black and white holes. Connecting an astrophysical black hole to a white 
hole, they lead to a natural mechanism of generating new universes. Under favorable conditions the 
new universes will also contain black holes which, in their turn, will give rise to another generation 
of universes. In this case the cosmological evolutionary tree will continue to grow to form the 
"hyperverse" . This scenario essentially differs from other known mechanisms, such as bounce, 
birth from "nothing", baby-universe scenario, etc. 
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I. INTRODUCTION 



Singularities in solutions of equations that are supposed to describe a physical phe- 
nomenon are usually associated with some kind of idealization. For example, in the case 
of hydrodynamics of ideal fluids one has multiple solutions in the problem of a body im- 
mersed in a flow a paradox resolved by taking into account a non-zero viscosity. As for 
gravitational singularities, they already appear in Newtonian gravity. A sphere of particles 
with initially equal radial velocities will collapse to a point with infinite density. However, 
this singularity disappears as soon as we introduce random perturbations to the initial ve- 



locities. In particular, effective cutoff o 
of dark matter halos may be the reason 



small-scale perturbations in N-body simulations 
2} why the singular Navarro-Frenk- White "cusp" 
profile |3| does not reproduce observations. 

In General Relativy (GR) the situation is somewhat different and the introduction of 
small perturbations in initial conditions does not prevent the formation of the black-hole 



singularity in quasi-spherical collapse [4j. Moreover, as it was recognized by a set of theo- 
rems [5J, the singularities in GR form under a wide range of initial conditions. A belief arose 
that it has to do with the idealization that is neglecting quantum-gravitational effects and, 
hence, corrections to the Einstein equations at energies of the order Mpi ~ 10 19 GeV. Indeed, 
a consistent treatment of quantum field theory in curved space-time makes it necessary to 
introduce corrections to the Einstein-Hilbert gravitational action (see for a review). Pos- 
sible consequences of such corrections for the problem of the black-hole singularity was first 
studied in [7]. Nevertheless, what is the exact form of the corrections is unclear until now, 
and the issue remains. 

The latter fact led to attempts of treating the problem phenomenologically. That is, all 
possible corrections are ascribed to an effective energy-momentum tensor (EMT) and the 
Einsteinian form of the equations if preserved. If we assume that the singularity (later on, 
we consider the Schwarzschild singularity) must somehow disappear than, at least near it, 
the effective EMT acquires an exotic form. One can go even further and put the principle 
of finiteness of curvature and/or density as a universal nature law |8j. Particularly, this 



principle led to a model [9j , in which the Schwarzschild geometry is attached to the de Sitter 
one in the vicinity of singularity where the squared Rimann tensor R a ^sR a/3lS achieves the 



Planckian value ~ Mpi- As the inflationary paradigm states [10(, the early universe had a 



quasi-de-Sitter symmetry. Hence, in the cosmological context the model gives a mechanism 
of generating a universe in the interior of a black hole. 

The idea that black holes can serve as generators of new universes was realized a long 



time ago [ll|]. It takes advantage of two remarkable features. First, a time-reversed solution 
of the Einstein equations is also a solution. Thus, in the interior of a white hole - the 
time-reversed counterpart of black hole - particles can only move away from the singularity, 
which is reminiscent of the geometry of an expanding universe. Second, the collapse into 
black hole ends with a singularity while our universe should have begun with a singularity 
(see 12| and refs. therein), and the two singularities may be, in fact, one. 

The phenomenological approach, however, does not permit to define the equation-of- 
state of the effective matter. It makes tempting to consider cases with a high degree of 
symmetry like that of the above-mentioned de Sitter space-time. In the meantime, there 
are indications that less symmetric structures of the singular region are possible. First, the 
black-hole interior exhibits the symmetry R x § 2 of the Kantowski-Sachs model. Second, in 
approaching to the singularity the Belinskii-Lifshitz-Khalatnikov regime with alternating 
Kasner eras is realized. One may ask the question: what is the advantage of this if the 
singularity persists? It does, but its structure may be less harmful in the sense that in some 
directions tidal forces remain finite. If the effective matter that forms the geometry moves 
only in these special directions, it can pass the singularity to form the geometry beyond it. 
In another context, this kind of situation ha ppe ns, for example, in the cosmological model 



with mixture of the Chaplygin gas and dust 



13| 



As applied to the black-hole singularity the described situation takes places for the class 
of the so-called integrable singularities proposed in [l^j]. The main assumption is that 
gravitational potentials in a Schwarzschild-like metric remain finite. Note that it does not 
imply the finiteness of the curvature tensor, though tidal forces on some geodesies are finite. 
This allows us to continue the effective-matter flow to the expanding white-hole region and, 
thus, realize a natural mechanism of generating new universes in the interior of a black hole 
without the extra assumption on the absence of the singularity. Since other known scenarios 
of cosmogenesis (bounce, birth from "nothing" , baby universes) use this extra assumption, 
our scenario is instrinsically different. 

In this paper we further study properties of space-time near the integrable singularity 
and generalize previously developed model of eternal black-white hole for the more realistic 



case of a universe born by an astrophysical black hole. The paper is organized as follow. 
In Sections [II] and III II we discuss the issue of continuation of Schwarzschild-like metrics 
through the integrable singularity and study spacetime structure near it. In Section HVl we 
classify geodesic behavior in its vicinity. Section |V] is a re'sume' of the model of black- white 
hole while in Section I VI I a model of an astrogenic (formed inside a collapsed astrophysical 
black hole) universe is presented. In Section IVHI we discuss possible physical nature of 
integrable singularities and in Section IVlIII - their role in the global picture of cosmogenesis. 
In Conclusion we briefly summarize basic properties of astrogenic universes. 

II. HOW TO CONTINUE THE SCHWARZSCHILD METRIC 

A black hole with positive external mass M > without rotation and charge is described 
in GR by the Schwarzschild metric in vacuum: 



where the variable r > is defined as the internal curvature radius of a closed homogeneous 
isotropic 2-space dY 2 = 7^ dy l dyi , which does not depend on y % (i = 1, 2), G is the gravita- 
tional constant. The F-space is invariant under the group of motions G3 (two translations 
that allow moving to an arbitrary point Y, and a turn around a point), and it can be reduced 
to the form dY 2 = r 2 dQ, where dQ = d6 2 + sin 2 ^^ 2 is the unit 2-sphere S 2 , y % = (0,(p). 
As the pole 9 = it is possible to take any of points S 2 , and the turn around the pole is 
described by the angle <p G [0, 2tt). 

The space X orthogonal to Y is given in the Eulerian gauge, in which one of the coordi- 
nates x coincides with r. Additionally, metric ([1]) is independent of the coordinate t 6 M 1 . 
Thus the black hole in vacuum has the group of motion G4 acting on the hypersurface 
C 3 = R 1 x §> 2 (three translations in C 3 + rotation in S 2 ). 

Topologically geometry (JT]) is a 4-cylinder with homogeneous 3-surface C 3 and radial 
coordinate r. Its reduction to positive values, r > 0, defines the Schwarzschild sector of a 
black (or white) hole. Attempts to extend this solution inevitably lead to regions occupied 
by matter. Therefore, the problem of analytical continuation of metric (CD) must be solved 
using more general GR metrics with matter, which we only restrict by requiring spherical 
symmetry G3 in S 2 . 




(1) 



r 
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Let us consider a class of such metrics, which in the orthogonal gauge (gu = 0) have the 
form 

do 2 

ds 2 = dX 2 - dY 2 = -MKdt 2 + ^L_-p (d9 2 + sin 2 # dip 2 ) , (2) 

where the real variable p~ l is not restricted by the sign and is defined as the F-space internal 
curvature scalar that does not depend on y l : 

= > ( 3 ) 

The 2-tensor Ricci is constructed from the metric 7^, which by spherical symmetry can 
always be reduced to the form dY 2 = p dQ. 

For general reference frames p,K.,Af are 4-scalars depending on coordinates of the X- 
space dX 2 = riudx I x J ' , where K is the kinetic term p and M is restricted by positive values: 

K = P -^— = -1 - 2$ , Af = N 2 > , (4) 
4p 

$ is the metric gravitational potential. The energy-momentum tensor has the form = 
diag(Tj , —pj_, —p±), where Tj = (e + p)uiu J — p8j. At p > functions e,p,p± and = 
(u 1 , 0, 0) describe the energy density, longitudinal and transversal tension, and 4-velocity of 
matter (uiu 1 = 1), respectively. Metric ([2]) reduces to solution (T1J in the part of the domain 
y/p = r > free of matter. Hence, by analogy with the vacuum case and by definition, the 
domain /C > will be later referred to as the T-region of space (j2j). There, in particular, the 



component T\ describes the pressure 
the sign of p. 

GR equations relate the metric and material scalars: 



15j, and the meaning of other components depends on 



$' = 2nGP , $ = 2nGTf , (5) 
W _ 2nG(E + P) 

N K 1 ' 

{pN^ P__{N^Y p_ (TP y (7) 

N +P± 2 _ AN*K ' ' W 



where the prime and dot denote the partial derivative with respect to p and t, respectively, 

By modeling the state of the effective matter, one can determine P and E from the Bianchi 
identities (J7J), and by integrating (JSJ) and (JB} over dx 1 from the external solution ([!]) into the 



future one can reconstruct the metric potentials. Similarly, in terms of the mass function 
m = m(x 7 ) equations (JSJ) read 

Gm ^ K dx J . . 

, s __, mj = iwpeiK Tf m , (9) 

where (j) = d/dx 1 and the absolute antisymmetric tensor in X has the form 

e u = | det (njj) ( ? ^ ) , e IJ = | det (n XJ ) |~ 1/2 ( ^ J ) . (10) 



III. SPACE-TIME NEAR THE INTEGRABLE SINGULARITY 

We are interested in solutions (j2J) generated by metric (JT|) that are geodesically complete. 
The sufficient condition for the metric to pass through p = is the fmiteness of potentials 
$ and iV and their derivatives. We call such models black-white (or black/white) holes 
with integrable singularity . They are described by the mean (averaged over the vacuum 
state (jl])) metrical 4-dimensional space without punctures, which provides a continuous ex- 
tension of affine parameters on world lines of test particles through the singular hypersurface 
p = 0. Using reference frames constructed on these particles we investigate the geometry of 
black-white holes outside the Schwarzschild sector. 

The effective matter generated by strong gravitational field near the singularity changes 
the space-time near p = and is the physical reason for the existence of integrable sin- 
gularity. The formation of matter in extremely strong gravitational fields is in accordance 
with the Le Chatelier principle: this is how nature reacts to a sharp increase in amplitudes 
of the metric potentials at p — > (which would continue to grow in the absence of mat- 
ter), thus, preventing them from divergence. Specific quantum-gravitational mechanisms 
of mutual transformations of material and gravitational degrees of freedom under extreme 
conditions need to be considered separately. Here, by developing papers [ijj], we postulate 
the continuity of the gravitational potentials ([2]) in presence of the effective matter. 

In the R-region of the space-time (j2J) with signature (+, — , — ,—), specified by the gener- 
ating metric ([T]), we have /C < 0, p = r 2 > 0. Here the hypersurface r = is degenerate and 
represents a time-like world line of the center of a spherically symmetric matter distribution 
of matter (for example, the center of a star). If there is no matter in the R-region, r > 2GM 
as follows from (TjQ), i.e. the hypersurface r = does not lie in the R-region. 



In the T-region we have K. > and the variable p can have any sign because the physics 
in this region is determined by non-linear quantum effects, and a priori there are no grounds 
to believe that the signature remains indefinite [if]]. Here the singular hypersurface p = 
splits 4-space in domains with different signatures (—, + , — ,—) for p > and (—, — , + ,+) 
for p < 0. The complete geometry depends on the distribution and properties of the effective 
matter near p = 0. 

We assume that the effective matter distribution maintains the symmetry of the generat- 
ing field. For example, the region that is evolutionary adjacent to the Schwarzschild metric 
(for example, T MV ^ for r < ro = const < 2GM and T^ v = at r > ro) conserves Killing 
t- vector and depends only on p. The region inside the star keeps the spherical symmetry and 
the field homogeneity of the star. We will use these constraints when constructing models 
in Sections EM 

Several properties of geometries ([2]) should be noted. 

• The continuity of the potential $ implies integrability of the function P(t, p) along 
lines t = const (see ([5])): 

$(t,p) = $o + 2ttG fpdp, $ = $(t,0). (11) 
Jo 

• The continuity of $ and N implies the integrability E(t,p) along lines t = const 
(see flSJ) 

N(t, p) = N exp (2irG J ^^dp^j , N = N{t, 0) . (12) 

• For /C = — 1 — 2$ = Kq = const > the structure of space-time at p — > has the 
form 

ds 2 = -dt 2 + dUdV - UVKl sin 2 fl dip 2 , (13) 
t = K jN dt, p = K 2 UV, d = ^- Q ln \^- 
The direction 9 can be taken along any meridional direction in § 2 counted from arbi- 
trarily chosen pole 6 = 0. 

Thus, the integrable singularity of a black-white hole includes hypersurfaces-horizons 
U = and V = lying in the T-region and intersecting along a space-like bifurcation line 
U = V = 0. They separate the T-region in cone sectors of the black (U < 0,V < 0) and 
white (U > 0, V > 0) holes and the static zone (UV < 0). The horizons have the cylindrical 



symmetry 



1 x § 2 with the space-like longitudinal axis i 6 I 1 and null geodesies in S 



2 



(t = const) 17| Photons propagating in these directions are kept gravitationally at p = 
by infinitely oscillating in Y within a finite interval of the affine parameter. Trajectories of 
other particles, including photons with projection in R 1 , intersect the singular hypersurface 
p = and escape in other metric domains. 

IV. GEOMETRIC MAPS OF AN OSCILLATING HOLE 

The structure of a 4-space can be described by 2- and 3-dimensional cross-sections which 
can be covered with a net of trajectories of test particles on them. It is convenient to choose 



coordinates such that the light trajectories have a slope of 7r/4 like in the flat world 18j. By 



spherical symmetry all geodesies in fl2]) can be represented by three groups of cross-sections: 

• X-planes - longitudinal cross-sections (9, <p) = const; 

• UV-pl&nes - transversal cross-sections (t, (p) = const; 

• tt/V-hypersurfaces - longitudinal-transversal cuts tp = const. 

The X-planes are filled with radial or longitudinal geodesies, the tUV-p\&nes contain parti- 
cles propagating in both longitudinal and transversal directions, and the UV-p\a.nes orthog- 
onal to the bifurcation line are filled with spiral geodesies which do not escape beyond the 
T-region (see Appendix) (Fig. [TJ. 

As an example, consider metric ©, independent of t, filled with a matter in the T-region 
r < To = const < 2GM. The matter is distributed in the t/V^-planes along bundles of 
straight lines (t, 9, <p>) = const intersecting at points U — V — 0. Integrating (jH]) along 
the time lines r = yfp from the black hole region (TjQ) into the future, we obtain continuous 
potentials 



Gm f f 

$ = , m = M — An p(r)r 2 dr = — 4rr / ppdr , (14) 

r Jr Jo 

which directly connect both holes. The black hole metric is in the sector r > 0, while 

the white hole is obtained from the extension of solution (|14|) into region r < under the 

finiteness condition pp = ^ Q /4:7iG = const at p = 0. The complete geometry can be restored 

by integrating from the bifurcation line points along all bundles of lines, including regions 

p < 0. 

8 




FIG. 1. The longitudinal (a) and transversal (b) cross-sections of an oscillating black- white hole. 
J + is the light infinity of the future for observers in the region r > 0, J~ is the light infinity of 
the past in the region r < 0, 1° are the space infinities of R-regions. 

For illustration, let us consider the model of an oscillating black-white hole with triggered 
matter distribution depending on p < p = r o ( see [141]' case B). In other terms, the 
transversal pressure increases in a jump: p± = Xq0{po — p). The longitudinal pressure, 
which is chosen to be vacuum-like (p = — e), is calculated using the Bianchi identities (J7j): 

p<p : P=-£=Ha„, * = l -H\ (p - 3p„) , (15) 

P>,o: P=-E = ^, » = -^, (16) 

47r | | |r| 

where Hf = 8irG\ /3 = GM/tq = const. By substituting ( fT5l) and (|T6l) into (j2D, we obtain 
the metric with matter (p < p ) 

4pK (1 + pf 

and the metric in vacuum (p = r 2 > p ) 

, 2 / 2GM\ , 2 dr 2 2 ^ 

rfs 2 = (l dt 2 2GM ~ ' ^ ^ 
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where i = K t, p = 4r 2 p(l + p)-\ K = K 2 - H 2 p, n = K*/H x = r (3 - r G /GMf/\ K 2 = 
— 1 — 2$o ; = ^(0) = —3GM/2rQ. The value of the potential of the integrable singularity 
is one and a half times lower than that at the matter boundary, and this ratio does not 
depend on the parameter r . 

Let us represent (JT7J) using the proper interval on the axis p: 

p = r 2 E [0, p ] : ds 2 = dr 2 - cos 2 (Fir) dt 2 - r\ sin 2 (Fir) dtt, (19) 

p < o : ds 2 = -dx 2 - cosh 2 (iJ 1 x) dP + r\ sinh 2 (if 1 x) rffi , (20) 

where r = —rism(Hxr) E [— r ,r ] and p = — r 2 sinh^ifxx) < 0. The space-time of static 
domains is asymptotically the anti-de-Sitter one: 

P«-r 2 : T^ = -A ^, W v = m\^ v . (21) 

For this reason the regions p < will be also referred to as anti-de-Sitter zones (ADS). 



V. THE SOURCE OF ETERNAL BLACK (WHITE) HOLE 

Eternal holes, in which matter is absent almost everywhere, can be obtained in the limit 
of small r — > in equations ffT5]) - ffT8]) : 

r E K 1 : e = -p = 2p ± = M ^ » ( 22 ) 

3M . . 

p < : p = -e = p ± = A = —j • 23 

87rrQ 

The ADS zones are strongly curved by the dense vacuum A (see (12"01 ): 

p « -r 2 : ds 2 = -dx 2 + \rle 2H ^ (-H 2 dt 2 + d9 2 + sin 2 ^<p 2 ) , (24) 

which produces a 5-like material source of geometry of the eternal black hole ffl8l) in the 
region r E R 1 . This source localized at |r| < r has the density ~ Ao and total mass M. 
The value ro can be calculated in quantum field theory (see an estimate in jl4]). In our 
classical treatment Tq is a free parameter of the problem. Note also that the relation between 
the longitudinal and transversal tensions in f l2"2"j) is not universal, it is related to the chosen 
condition of the model p = — e. 

Solution ( 122|) -( )24|) shows that the polarized vacuum in static ADS zones is the source of 
eternal black holes. The gravitational mass of the effective matter of each of the ADS-zones 
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generates, via its light hypersurface of the future p = 0, the causally connected Schwarzschild 
metric of a black-white hole in vacuum, which extends in time up to the next integrable 
singularity p = 0, where the process repeats (see Fig. [I]). The full geometry is invariant with 
respect to the inversion r — > — r of the spherical coordinate system relative to the bifurcation 
lines. In this sense phase transitions between gravitational and material degrees of freedom 
in this model are reversible. 

In Section [VI] we construct one more example of reversible geometry, in which the black- 
hole source is a parent star that collapses from the R-region with the formation of a singu- 
larity in the T-region. In the course of evolution the space-time passes the stage of effective 
matter and transforms into a white hole with the metric of a homogeneous cosmological 
model. 

VI. ASTROPHYSICAL BLACK- WHITE HOLE 

Let us consider the model of a black-white hole with integrable singularity assuming 
that the black hole was formed during the collapse of a star from the parent universe. The 
star is modeled as a homogeneous sphere with radius r(T,R = 1) = a(T), which was at 
rest in a flat space-time and had the mass M (a ^> 2GM at T — > — oo) and then started 
collapsing due to self-gravitation, with initial pressure being negligible. Here T and R are 
the Lagrange coordinates co-moving with the star material, and the proper time T e M 1 
and radial markers of spherical shells R > are normalized such that R = 1 at the star's 
surface. 

The field symmetry is determined by the initial conditions of the problem. Inside the 
star (R < 1) the symmetry IR 1 x E 3 with 6-parametric group G$ on E 3 and potentials 
a, H, e,p are functions of T. Outside the star (R > 1) the group of motion G A on C 3 and all 
potentials depend on r. To avoid crossing of matter flows we assume that the longitudinal 
tension, which is produced outside the star at large curvatures, is vacuum-like. Then at 
R > 1, N = 1 and the energy- momentum tensor is invariant with respect to motions in X: 
Tjf = —diag(p,p,p±,pj_). Inside R < 1 the matter is Pascalian, = (e + p)u u u >l — p5^ 
(where 4- velocity u v = T u ), and its state is calculated from boundary conditions. The 
homogeneity of the star suggests the continuity of $ and p at the boundary R = 1 (r = a), 
whereas the density and transversal pressure can change discontinuously. 
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Inside the star the metric has the form 12j 



ds 2 = dT 2 - a 2 (dR 2 + R 2 dQ) = - r 2 dQ , (25) 

v / 1 + 2$ 

where the Euler and Lagrange coordinates are related as 

f dT R 2 

R<1: r = aR, t = --—- — , (26) 



a?H 2 

and the Hubble function H = da/adT can be found from the Friedmann equations: 
tt9 8irG 2Gm 2$ de , 

H =—t = — = -^> ^ + 3H(e + p) = 0. (27) 

Initially the pressure is absent and the star collapses freely (H < 0, a 3 H 2 = 2GM). The 
Newtonian potential in this limit is $a? = GM(R 2 — 3)/2a. The internal tension in the star 
appears at a < r < 2GM, and it can be calculated by sewing it to the effective matter at 
the boundary R = 1. 

Outside the star the metric has the form 

dv 2 

ds 2 = (1 + 2<5>)dt 2 - i + 2$ - r 2 dtt = dT 2 + 2$dR 2 - r 2 dVt , (28) 

where the Lagrange reference frame co-moves with shells of free dust particles following the 
collapse of the star surface: 



f dr f V^2$dr . . 

R>1: R-T= ■ , t = T- — , 29 

J V=2$ J 1 + 2$ K J 

Gm 1 f 

$ = = — H 2 r 2 , m = M-4n / p(r)r 2 dr. (30) 

r 2 J ro 

The source of the Schwarzschild metric is the mass M of the star with the continuous 

potential $ = —GM/r at the star boundary. The discontinuity of functions N, t and gun is 

due to the density jump at R = 1. The effective matter outside the star emerges at r < r 

and preserves the symmetry of the parent metric. The continuity of $ suggests m(0) = 

and leads to formulas of Section IIV1 

By extending the reference frame of free particles over the entire cone p = r 2 > 0, we 

continue metric from the black hole region (r > 0) into the white hole region (r < 0) 

bypassing the body of the star itself (R > 1) (see (I14p - ([I5]) ). Closing the external metric, 

we then restore the complete solution on the entire manifold R > 0, by sewing p at the star 
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FIG. 2. The Penrose diagram of an astrophysical black-white hole. The hell gray region is the 
body of the star, the dark gray region is the effective matter, the dashed line separates the R- and 
T-regions inside the star, the line with arrows is the contour t = const. 

boundary (Fig. [2]). As a result, everywhere in the effective matter zone we obtain P = 2Ao/3 
and 

\r\<r R: r = -VSroRsm^f ) , $ = ^%R 2 cos 2 (HiT) , (31) 

where R = min(i?, 1), T = T + (1 — R) ■ 6(R — 1). The state of effective matter inside the 
star is described by the equation e + 3p = 2X . 
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Freely falling shells of the star intersect simultaneously in one point r = separating the 
time-like and space-like parts of the axis r = (see Fig. [2]). At this bifurcation point the 
equation of state of the effective matter is p = — e/3 and the motion is rectilinear: a oc T. 
Gravitation is 'switched off' due to the zero mass at r = 0, which provides finite amplitude 
of the tidal forces and the continuation of world lines of the matter into future. 



VII. PHYSICAL NATURE OF INTEGRABLE SINGULARITIES 

The concept of integrable singularities allowed us to advance in solving the cosmogenesis 
problem using the new class of black-white hole models in GR. Are these solutions only 
research tools or they really exist? If yes, then why and how do they form? These questions 
require separate studies. Here we discuss some ideas on the physical nature of integrable 
singularities. 

Reasons for the formation of singular structures with finite gravitational potential are 
related to the redistribution of material and space-time degrees of freedom in strong grav- 
itational fields of singularities themselves. (At high intensities and small scales quantum 
effects apparently play the key role). That is why we speak of the effective matter, which 
includes both material and gravitational degrees of freedom that remain after the metric 
averaging over the quantum state. However, we keep notions of energy-momentum and laws 
of motion in the form of the Bianchi identities. This allows us to describe the back reaction 
using the Einstein equations in which the left side is defined using the mean metric g^, and 
the right side contains the energy-momentum tensor of the effective matter including all 
polarization, gravitationally modified and other terms of quantum theory that need to be 
calculated. 

Let us illustrate the method by considering the behavior of one degree of freedom in the 
given physical symmetry of the collapsing star. Let the physical variable be described by 
a massless field p minimally coupled to metric ( 12"5|) . The second quanization leads to the 



following equation for amplitudes of Fourier harmonics p/~ — Vk{ri)/a 12] : 



v»+(k 2 -U)v = 0, (32) 
U = U(rj) = — = a 2 H 2 + {aH)' = -2$ + -^== , (33) 
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here the prime denotes the derivative with respect to the conformal time rj = J dT/a, k is 
the wave number, $ is the gravitational potential on the stellar surface. 

As long as the pressure in the star is low, the tidal potential increases with decreasing 
radius: U = GM/a = a 2 H 2 /2. At large a the state of the field is vacuum-like: it oscillates 
near an equilibrium point in the adiabatic zone U < k 2 . During the collapse the field enters 
the parametric zone, and its amplitude increases unlimitedly as U > k 2 increases: 



aV2k k 6 ' 1 



At this stage, the pressure cannot be neglected any more, vacuum is polarized and, in its 
turn, affects the metric, which changes the gravitational potential and the rate of collapse 
near a = (for the solar mass and the number of degrees of freedom ~ 100 this occurs at 
H ~ 0, 1 x M P (\nU y 1/2 ~ 10 17 GeV). 

One can suppose that the back reaction, described by the Friedmann equations, restruc- 
tures the solution in such a way that tidal potential fl33|) stops increasing and saturates: 
U(r — > 0) = Uq = — 2$ — const. This precludes the ultraviolet catastrophe and saves 
oscillators from being destroyed, because high-frequency modes with k 2 > Uq remain in the 
adiabatic zone and are not polarized. Thus the notion of the Heisenberg state vector that 
stems from the Minkowski vacuum in the R-region stays valid, and the main requirement of 
the metric theory on the finiteness of the gravitational potential on the singular hypersurface 
is secured. We call such a structure an integrable singularity. 

This example illustrates the difference of our theory from models with bounces (see, for 



example, [llf|), which have no singularity, a > 0, and at the bounce H = $ = 0. The last 
requirement, in our opinion, is redundant. In our models the potential $ achieves extremum 
at the singularity r = 0, which weakens the singularity but does not fully eliminate it. For 
small a oc T the situation resembles the explosive model of non-gravitating particles by 
Milne. At the moment of crossing, particles move by inertia and do not feel the central mass 
attraction. Gravitation at this moment is switched off (m(0) = 0), therefore the space does 
not bend, although the energy density diverges. 
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VIII. THE ARROW OF TIME AND NATURAL SELECTION 



Gravitationally tidal models of black-white holes containing expanding flows of matter 
elucidate many fundamental unanswered questions in the modern physics. One of them is 
the causality principle, which puts the causal relationships in correspondence with the arrow 
of time. In the context of our work we can discuss the origin of the cosmological arrow of 
time, which we consider as the orientation of the future light cone in the direction of the 
volume expansion of a large-scale flow of matter. 

As is well known, dynamical equations describing microscopic processes are invariant 
relative to the change of the sign of time. However, the local dynamics has to be completed 
with an external arrow of time, since the invariance relative to the direction of time is lost 
in the limit transition to the global geometry. 

Geodesically full geometries with integrable singularities suggest the origin of the arrow 
of time. They include different space-time domains separated by event horizons r = and 
r = 2GM: non-stationary regions of black and white holes (lying between the Schwarzschild 
and singular horizons) and alternating static R- zones (connecting the Schwarzschild horizons 
of white and black holes), and ADS-zones (connecting singular horizons of black and white 
holes). Here all possibilities are realized. Each collapsing and anti-collapsing (cosmological) 
region has its own time. The static regions are time-independent but are not spatially 
homogeneous. 

We deal with a unique geometry that is split into sectors with time and space parity. 
When crossing any of the horizons the sense of the coordinate r, which the metric depends 



on, changes 19]. In some domains r is the time coordinate (and then we obtain black holes 
and/or cosmological models), in other domains it is the space coordinate (static zones). The 
complete geometry here can remain invariant relative to the change r — > — r. 

Therefore, the origin of the cosmological arrow of time is related to the initial conditions. 
We (observers) belong to a cosmological flow of matter and live in its proper time which 
started at the moment r = some 14 billion years ago. The time coordinate can be 
extrapolated into the past, to the pre-cosmo logical epoch, but there it described the time 
on the T-zone of the parent black hole. Furthermore, this coordinate represented the radial 
coordinate of an asymptotically flat space of the parent universe, in which the star had lived 
before it collapsed to the black hole. The integrable singularity r = that emerged during 
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the collapse 'kindled' our Universe, and after several billion years the non-linear large-scale 
structure began developing to initiate the process of star formation. As a result, new black 
holes appeared; they can be entrances to new universes. 

This process of evolution of the multisheet space-time resembles the growth of a tree (a 
genealogical tree, so to speak). Such a tree can flourish, or can wither if no new black holes 
are formed in the daughter universes. The critical situation appears when the development 
conditions do not provide the production of seed density fluctuations to form gravitationally 
bound matter clumps and their collapse to black holes. But another scenario can be realized: 
one collapse under favorable conditions, which realize inflationary parameters and phase 
transitions in a white hole, can lead to flourishing of a whole tree with non-decaying chains 
of new universes. Because of these processes the cosmological natural selection works 20]: 
only universes where black holes can be formed survive and develop, and this is possible for 
a certain set of parameters, world constants, etc. 

This concept of multisheet hyperverse is based on the gravitational instability processes 
which resemble oscillating tides. Anti-collapsing space-time regions (white holes) stem from 
collapsing black holes, and, inversely, an expanding quasi-homogeneous flow of matter of 
a white hole desintegrates into clumps collapsing into black holes. The former process 
is related to the r = horizon, while the latter - to the r = 2GM horizon. At both 
horizons the gravitational potential is relativistic, and quantum-gravitational processes of 
vacuum polarization and pair creation should be taken into account. However, while on 
the Schwarzschild horizon these effects are suppressed by the mass parameter (the Hawking 
evaporation), at r = they dominate and form structures of integrable singularities. 



IX. CONCLUSIONS 

In the framework of our concept of integrable singularities a new class of black- white holes 
in GR is obtained, which can be the key to solve the cosmogenesis problem. The integrable 



singularity r = can be compared with a classical cusp in which 



;he energy density or 

the longitudinal tension of matter diverge, but the mass is zero [211 ]. m(0) = 0, and the 
gravitational potential $o is limited. Because of this property the tidal forces are finite 
and the geodesies freely extend from a black hole into a white hole, where the geometry is 
equivalent to that of an expanding cosmological model. 
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The mass of matter in the multisheet hyperverse can be arbitrarily large, since it is 
compensated by the negative gravitational bound energy. Thus the total energy of holes, 
measured in static zones, is constant in time. The integrable singularities is reminiscent 
of machines for reprocessing gravitational degrees of freedom into material ones, however 
quantitative characteristics of this process can be determined only in the self-consistent 
quantum theory. 
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Appendix: Motion of test particles in a black-white hole 

Consider the motion of test particles in metric (jSJ), which is independent of coordinate t. 
Let k^ = dx^/dX and A be the tangent vector, and the affine parameter along the trajectory 
x ti _ x v>(\y From the geodesic equation fc M /c„ ;At = we obtain 

k 

k = k t = const , k± = kg H ^— = const , (A.l) 

sin 9 

where k v = const is the azimuthal angular momentum. Unlike invariants of motion for 
the longitudinal and transversal momenta, the value of the azimuthal number depends on 
the orientation of the polar coordinate system. By adjusting the pole 6 = with one of 
the points on the trajectory, we have ip = const and k^ = everywhere on the world line 
of a particle. In the projection onto a 2-sphere the particle moves along a meridian with 
monotonic increase of the angle 9 within the interval 2tc by the number of turns in § 2 : 




(A.2) 



The equation for X(p) follows from the normalization integral k^ = n = const, where 
n = for the null and n = 1 for the time-like geodesies, respectively. Inside the cone 
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FIG. 3. The potential 4>{r) of light geodesies in a black-white hole. 

p = r 2 > of a black-white hole equation ( 1A.2I) has the form: 
/ dr\ 2 k 2 k 2 

[^) +0 = ^-n, 0=-f (l + 2$) + 2n$, (A.3) 

where the potential <fi = <ft(r) tends to zero at r — > oo (Fig. [3]). 

The longitudinal light (n = k± = 0) geodesies propagate directly from the black hole 
into the white hole with continuous affine time A = —k~ l J Ndr. The spiral light geodesies 
on the cross-cuts (n = k = 0) lie in the T- region of the black or white hole, K. = K 2 > 0, 
coming from ADS-zones and going back to ADS-zones: 9 = —k± J dX/p = ± J dp/(2pK). 
Of special interest are photons with n = k = k± = living on the horizons p = 0: 

g = ± ln|A-A | ^ p = T 2K (\-\ )k ± ^0, (A.4) 

where Ao is the value of the affine parameter on the bifurcation line, 

Trajectories of photons with the impact parameter k = /cj_(3a/3GM/c) _1 < 1 connect 
both R- and T-regions. Photons with k, = 1 are turning in R-regions at radius r = 3GM. 
For k, > 1 there are photons of two kinds: in the R-region with r > 3GM, and in the zone 
r < 3GM which unites the T-region and the inner part of the R-region adjacent to it. 
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